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ABSTRACT 


The  eclipsing  binary  system  UV  Piscium,  which 
is  a  short-period  RS  Canum  Venaticorum-type  binary  system, 
was  observed  at  Mt.  Laguna  Observatory  during  1980  using 
the  Stromgren  4-color  filter  system,  and  broad  and  narrow 
Calcium  II  K-line  .filters.  A  new  ephemeris  was  calculated 
using  only  photoelectric  times  of  minima.  Four  sets 
of  orbital  solutions  were  obtained  using  the  Nelson-Davis- 
Etzel  model  giving  reliable  geometric  elements  for  the 
system.  The  orbital  separation,  radii  and  masses  of 
the  stellar  components  are  given  in  absolute  units  based 
on  Popper's  spectroscopic  masses.  The  light  curves 
exhibit  small  distortions  in  the  maxima  light  levels 
which  vary  with  phase.  The  system  is  composed  of  two 
well -detached,  main  sequence  stars  with  the  primary  a 
G2  star  and  the  secondary  an  early  K  star. 


ACKNOWLEDGEMENTS 


I  am  highly  appreciative  of  the  many  people  who 
aided  me  in  the  completion  of  this  study.  Foremost  among 
them  is  Dr.  Burt  Nelson  for  his  original  suggestions 
leading  to  this  thesis  topic  and  his  helpful  advice 
throughout  the  course  of  the  work  on  it.  Appreciation 
is  also  extended  to  Dr.  Fred  Talbert  for  his  many  helpful 
suggestions  and  advice,  and  to  Dr.  Don  Rehfuss  for  serving 
on  my  thesis  committee. 

Special  thanks  is  extended  to  all  the  assistants, 
especially  to  JoAnn  Eder,  who  aided  me  in  the  completion 
of  the  observations  at  Mt.  Laguna.  Without  their  help 
the  observations  for  this  project  would  not  have  been 
obtained. 

Finally,  my  greatest  thanks  goes  to  my  wife, 

Irma,  and  my  family,  whose  support  and  understanding 
throughout  the  many  nights  at  the  observatory  and  the 
computer  facilities  is  lovingly  appreciated. 


1 


TABLE  OF  CONTENTS 


Page 


ACKNOWLEDGEMENTS . iii 

LIST  OF  TABLES .  v 

LIST  OF  FIGURES .  vi 

Chapter 

1.  INTRODUCTION  .  1 

RS  Canum  Venaticorum  Binaries  .  1 

Historical  Background  on  UV  Piscium  ....  4 

2.  OBSERVATIONAL  PROGRAM  .  9 

3.  EPHEMERIS  CORRECTION .  18 

4.  LIGHT  CURVES  AND  THEIR  SOLUTION  .  26 

Light  Curves .  26 

Solution  of  the  Light  Curve .  39 

5.  CONCLUSION . .  .  46 

REFERENCES .  50 

APPENDICES .  53 

A.  Photoelectric  Data . 54 

B.  Timings. of  Primary  Minimum .  74 

ABSTRACT . 76 


MMtoi m 


LIST  OF  TABLES 

Table  Page 

1.  Catalog  Data .  12 

2.  Comparison  Star  Transformation  Values  ....  17 

3.  Photometric  Minima  of  UV  Piscium .  20 

4.  Model  Parameters .  41 

5.  Solution  Values  and  Standard  Errors  .  43 

6.  Mean  Values  and  Standard  Deviations .  44 

7.  Absolute  Parameters  .  46 


7.  &(b-y)  Curve .  33 

8.  A(u-b)  Curve .  34 

9.  Calcium  II  K-Line  for  UV  Piscium .  35 

10.  Calcium  II  K-Line  for  Comparison .  36 


Chapter  1 


INTRODUCTION 

This  study  was  undertaken  in  order  to  add  to 
the  growing  accumulation  of  observational  data  on  the 
RS  Canum  Venaticorum  type  eclipsing  binary  systems. 

This  class  of  binaries  has  only  in  the  last  decade  become 
recognized  as  a  distinct  and  separate  set  of  binary 
systems  with  its  own  distinct  set  of  physical  properties. 
Many  of  the  systems  identified  as  such  have  very  sparse 
accumulations  of  data,  particularly  in  the  photometric 
area  of  observation. 

RS  Canum  Venaticorum  Binaries 

The  RS  Canum  Venaticorum  binaries  have  been  super 
ficialiy  referred  to  as  a  subset  of  Algol-like  eclipsing 
binaries  but  with  some  unusual  properties  which  differen 
tiate  them  from  normal  Algols  (Morgan  and  Eggleton  1979'. 
They  were  specifically  defined  in  a  review  paper  by  Hall 
(1976)  based  on  their  observed  properties  and  physical 
characteristics.  To  be  considered  a  member,  a  binary 
system  must  have  an  orbital  period  between  one  day  and 
two  weeks,  a  hotter  component  of  spectral  type  F  or  G 
with  luminosity  class  V  or  IV,  and  H  and  K  emission  in 
the  spectrum  outside  eclipse.  Hall  listed  24  systems 


which  met  these  criteria. 

There  are  at  least  13  properties  characteristic 
of  the  RS  Canum  Venaticorum  binaries,  although  no  one 
system  necessarily  exhibits  all  of  them.  (1)  The  H  and 
K  emission  is  from  the  cooler  star  (or  both  stars). 

(2)  The  cooler  star  is  around  a  spectral-luminosity  class 
of  KO  IV.  (3)  Ho  emission  is  seen  outside  of  eclipse. 

(4)  A  wave-like  distortion  is  present  in  the  light  curve 
outside  of  eclipse.  (5)  This  wave  migrates  toward 
decreasing  phase.  (6)  Additional  irregular  variations 
not  due  to  the  wave  are  present  in  the  light  curve. 

(7)  One  or  both  components  has  an  ultraviolet  excess. 

(8)  One  or  both  components  has  an  infrared  excess.  (9) 
Strong  centimeter  radio  emission  has  been  observed  in 
some  systems.  (10)  Very  large  and  irregular  changes 
in  the  orbital  period  have  been  observed  which  cannot 

be  explained  by  apsidal  motion  or  a  third  orbiting  body. 
(11)  Some  period  variations  can  be  correlated  with  the 
migrating  wave.  (12)  Most  systems  have  a  mass  ratio 
near  unity.  (13)  Both  components  are  detached  from  their 
respective  Roche  lobes. 

The  most  conspicuous  of  these  properties  and 
the  one  which  has  attracted  the  most  attention  is  the 
distortion  wave  in  the  light  curve.  The  wave  is  typically 
a  persistent,  nearly  sinusoidal  wave  which  distorts  the 
light  curve  both  inside  and  outside  of  eclipse  and  renders 
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the  two  maxima  unequal  in  brightness.  The  wave  has  a 
period  nearly  equal  to  the  system  orbital  period,  but 
migrates  slowly  toward  decreasing  orbital  phase.  Oliver 
(1974)  was  the  first  to  point  out  that  these  waves  and 
their  migration  seem  to  be  common  to  the  RS  Canum  Venati- 
corum  binaries  in  general.  The  migration  periods  have 
been  found  to  range  from  5  to  75  years,  and  the  wave 
amplitude  also  fluctuates  in  cycles  with  time  scales 
of  decades  (Hall  1976).  The  wave  can  be  attributed  to 
the  cooler  star  (Biermann  and  Hall  1976),  and  the 
migration  of  it  can  account  for  the  variable  depth  of 
the  primary  minimum  and  the  variable  displacement  of 
the  secondary  minimum. 

Several  different  explanations  have  been  advanced 
to  explain  these  distortions.  One  was  suggested  by 
Catalano  and  Rodono  (1974)  in  which  the  precession  of 
a  ring  or  disk  of  circumstellar  matter  tilted  with  respect 
to  the  orbital  plane  produces  the  distortion  wave.  A 
second  was  expanded  on  by  Popper  (1977)  which  would 
account  for  the  distortion  by  suggesting  one  of  the 
stellar  components  is  pulsating. 

The  third  explanation,  the  star  spot  hypothesis, 
was  actually  the  first  one  proposed  (Kron  1947)  and  has 
regained  prominence  recently  (Hall  1972;  Eaton  and  Hall 
1979).  This  concept  has  an  advantage  in  that  it  organizes 
the  stellar  surface  activity  analogously  to  that  observed 


on  the  Sun  but  on  a  greater  scale.  The  concept  of  an 
active  chromosphere  also  provides  a  possible  explanation 
of  many  of  the  other  observed  properties  of  these  systems 
including  the  Ho  emission,  H  and  K  emission,  and  period 
changes . 

In  addition  to  the  24  systems  originally  listed 
by  Hall  (1976)  as  RS  Canum  Venaticorum  type  binaries, 
to  which  additional  systems  have  since  been  added,  Hall 
also  listed  two  subgroups  which  would  be  classified  as 
RS  Canum  Venaticorum  type  systems  except  for  their 
periods,  which  fall  outside  of  the  defined  range.  These 
are  the  long-period  and  short-period  groups.  These  sys¬ 
tems  exhibit  many  or  all  of  the  other  characteristics 
common  to  the  actual  RS  Canum  Venaticorum  binaries. 

Hall  (1976)  specifically  defines  the  short-period  group 
as  non-contact  binaries  with  periods  less  than  one  day 
in  which  the  hotter  component  is  of  spectral-luminosity 
class  F-G  V-IV,  and  H  and  K  emission  is  displayed  in 
one  or  both  components.  Hall  listed  six  systems  in  this 
group  of  which  UV  Piscium  was  one. 

Historical  Background  on  UV  Piscium 

UV  Piscium  is  a  ninth  magnitude  eclipsing  binary 


star  system  with  a  period  of  0.86  days.  This  system 
was  selected  as  having  a  period  and  observing  season 
which  would  allow  complete  coverage  of  the  light  curve 


with  the  41-cm  Boiler  and  Chivens  telescope  operated 
by  San  Diego  State  University  at  Mount  Laguna  Observatory 
during  the  Fall  of  1980. 

The  variability  of  UV  Piscium  was  first  discovered 
in  1957  at  the  Remeis-Bamberg  Observatory  where  it  was 
designated  BV  149.  The  first  light  curve  was  published 
by  Huth  (1959)  based  on  348  photographic  plates.  Huth 
gave  a  maximum  light  magnitude  of  9.6  with  a  range  of 
0.9  and  0.2  magnitudes  for  the  primary  and  secondary 
minima  respectively ,  but  no  orbital  solution  was  attempted. 
The  period  derived  from  his  minima  was  0.861046  days. 

Subsequent  light  curves  were  investigated  by 
Carr  (1967)  and  Oliver  (1974).  Carr's  light  curve 
observed  in  1966  was  the  first  to  be  obtained  by  photo¬ 
electric  means.  He  concluded  the  three  light  curves 
obtained  using  UBV  filters  resembled  an  Algol  type  system. 
The  light  curves  contained  an  asymmetry  in  the  light 
levels  of  0.05  magnitudes  in  the  two  maxima  at  phases 
0.25  and  0.75.  Carr's  solution  of  the  light  curves  indi¬ 
cated  that  the  eclipses  were  complete,  that  the  primary 
eclipse  was  a  transit  and  that  the  secondary  was  an  occul- 
tation.  His  solution  was  obtained  by  using  the  procedures 
of  Russell  and  Merrill  (1952). 

Oliver's  light  curve  of  UV  Piscium  was  obtained 
in  1968  and  1969  in  conjunction  with  observation  of  twenty 
eclipsing  binary  systems  which  exhibit  Calcium  II 
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emission.  Unfortunately  his  light  curve  was  incomplete 
with  no  observations  of  the  maxima  between  phases  0.15 
and  0.40.  In  addition,  his  solution  was  obtained  using 
the  nomographs  due  to  Merrill  (1950)  and  was  at  odds 
with  that  of  Carr.  Oliver  concluded  that  the  eclipses 
were  partial  and  that  the  primary  star  had  the  smaller 
radius.  This  conclusion  was  reached  despite  the  fact 
that  a  visual  inspection  of  his  secondary  eclipse  obser¬ 
vations  indicates  a  flat  bottom  which  would  suggest  a 
total  eclipse  due  to  an  occultation. 

Sadik  (1979)  published  a  comprehensive  analysis 
of  UV  Piscium  based  on  UBV  photometry  observed  in  1977. 

He  was  able  to  obtain  absolute  elements  for  the  system 
using  Popper’s  spectroscopic  elements  (1976).  His  solu¬ 
tion  was  consistent  with  Carr's  in  that  the  primary 
eclipse  was  due  to  a  transit  and  the  secondary  due  to 
an  occultation.  The  primary  was  the  larger  of  the  two 
stars.  Sadik  assumed  the  primary  to  be  a  G2  V  and  the 
secondary  a  K0  IV. 

The  assumption  of  spectral  types  used  by  Sadik 
was  apparently  based  on  those  suggested  by  Oliver  (1974) 
which  were  derived  from  a  color  analysis  of  the  system 
and  the  relative  sizes  of  the  two  components  obtained 
by  Oliver  in  his  study.  The  assumption  of  a  subgiant 
class  for  the  secondary  when  this  star  has  a  smaller  radius 
than  the  primary  is  in  conflict  with  the  main  sequence 
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classification  of  the  primary.  According  to  the  stellar 
radii  listed  in  Landolt-Bornstein  (Voigt  1965),  a  KO 
IV  star  has  a  typical  radius  of  3.0  solar  radii,  and 
a  KO  V  star  has  a  radius  of  only  0.85  solar  radii.  Sadik 
gives  the  radius  of  the  secondary  star  as  0.929  solar 
radii,  which  would  indicate  it  should  be  considered  a 
main  sequence  star.  Popper  (1976)  classified  the  primary 
as  a  G2  and  gives  a  (B-V)  value  for  the  secondary  of 
+0.91  without  any  luminosity  classification.  This  (B-V) 
value  is  that  derived  by  Oliver. 

Sadik  also  detected  an  asymmetry  in  the  light 
curve.  He  concluded  that  the  asymmetry  could  be  due 
to  either  a  dark  spot  or  a  hot  spot  on  the  secondary. 
However,  a  hot  spot  was  determined  to  more  adequately 
explain  the  color  dependency  of  the  light  curve  asymmetry. 

Based  on  these  previous  studies  of  UV  Piscium  and 
the  current  work  on  RS  Canum  Venaticorum  systems,  an 
observational  program  was  devised  using  the  4-color  ubvy 
filter  system.  The  object  of  the  program  was  to  obtain 
three  complete  light  curves  in  each  filter  for  the  months 
of  September,  October  and  November,  1980.  From  these 
observations  a  new  solution  to  the  light  curve  could  be 
obtained,  as  well  as  an  examination  for  short  term  asym¬ 
metrical  variations  in  the  light  curves  from  month  to 
month.  Better  information  for  photometric  classification 
of  the  luminosity  classes  of  the  stars  was  also  sought. 
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Observations  were  also  planned  using  broad  and  narrow 
Calcium  II  K-line  filters  in  an  effort  to  detect  any 
phase  related  emission  effects. 


Chapter  2 


OBSERVATIONAL  PROGRAM 

The  photometric  observations  of  UV  Piscium  were 
obtained  on  12  nights  between  6  September  and  14  December, 
1980  U.T.  Two  additional  nights  of  ubvy  standard  star 
observations  were  also  obtained  in  this  same  period. 

These  additional  observations  were  used  to  determine 
transformation  coefficients  in  order  to  tie  the  comparison 
star  into  the  standard  system. 

All  of  the  observations  were  obtained  with  the 
41-cm  Boiler  and  Chivens  Cassegrain  reflector  telescope 
at  the  Mount  Laguna  Observatory  operated  by  the  San  Diego 
State  University  Department  of  Astronomy.  The  photometric 
system  consisted  of  a  single  channel  photometer,  incor¬ 
porating  a  dry  ice  cooled  1P21  photomultiplier.  A  37 
arc  second  diameter  diaphragm  was  used  on  all  observations 
except  for  the  first  night  when  a  72  arc  second  diameter 
diaphragm  was  used.  This  night  was  the  first  of  the 
two  standard  star  observation  nights.  The  filter  system 
consisted  of  the  4-color  uvby  interference  filter  system 
of  Stromgren  (1966)  and  a  broad  and  narrow  band  filter 
set  centered  on  the  Calcium  II  K-line. 

The  photomultiplier  signal  was  monitored  by  a 
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Weibrecht  charge  integration  amplifier.  Integration 
times  of  10  seconds  were  used  on  the  standard  stars, 
but  times  of  40  to  50  seconds  depending  on  the  filter 
used  were  required  for  the  program  and  comparison  stars 
in  order  to  insure  good  photon  statistics  and  to  keep 
the  gain  settings  in  a  range  where  linearity  between 
settings  was  maintained.  Data  recording  was  obtained 
with  an  automated  digital  system  consisting  of  a  digital 
voltmeter  accurate  to  four  significant  digits,  a  digital 
clock  set  by  WWV  radio  time  signals  to  the  nearest  second, 
an  internal  timer,  various  code  switches,  a  parallel 
to  serial  data  converter  and  a  teletype  equipped  with 
a  paper  tape  punch  unit. 

Following  completion  of  the  observing  program, 
the  filters  were  measured  on  a  Perkins-Elmer  double  beam 
spectrophotometer  (Coleman  124)  belonging  to  the  San 
Diego  State  University  Department  of  Natural  Science. 

The  mean  wave-lengths  of  the  u,  v,  b  and  y  filters  were 
found  to  be  3440,  4120,  4712,  and  5512  &  with  bandwidths 
of  387,  177,  172  and  232  X  respectively.  These  values 
proved  tp  be  a  relatively  close  match  to  the  filter  system 
used  by  Crawford  and  Barnes  (1970)  to  define  the  standard 
system  for  uvby  photometry.  Their  mean  wavelengths  were 
3500,  4100,  4700  and  5480  X  with  bandwidths  of  380,  200, 
200  and  200  X  respectively. 

The  broad  K-line  filter  had  a  mean  wavelength 
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of  3031  X  and  a  bandwidth  of  87  X.  The  narrow  band  In¬ 
line  filter  was  found  to  have  a  mean  wavelength  at  3940 
X,  but  it  also  had  a  secondary  transmission  peak  of  nearly 
equal  strength  at  4000  X.  The  narrow  band  filter  also 
had  a  bandwidth  of  96  X,  which  exceeds  that  of  the  broad 
band  filter.  Hence,  any  results  obtained  from  the  K- 
line  filters  were  held  to  be  questionable  at  best.  In 
addition,  the  K-line  filters  were  only  used  during 
September  and  October,  1980,  after  which  they  were  not 
used  in  order  to  obtain  greater  coverage  of  the  light 
curve  with  the  4-color  uvby  filters. 

The  technique  of  differential  photometry  was 
employed  for  the  program  star.  The  observations  were 
normally  made  repeatedly  in  the  order  sky,  variable, 
comparison  with  one  deflection  in  each  of  the  six  filters 
in  the  order  y,  b,  v,  u,  broad  K-line,  narrow  K-line. 
Occasionally  an  additional  reading  for  a  filter,  particu¬ 
larly  the  u  filter,  was  obtained  to  monitor  for  any  large 
scale  atmospheric  variations  such  as  cirrus  clouds. 

Care  was  taken  to  insure  the  program  star  was  bracketed 
between  sky  and  comparison  observations.  Observations 
of  the  check  star  were  interspersed  with  the  comparison 
star  each  night  in  order  to  monitor  for  any  variability 
of  the  comparison  star.  No  variability  was  detected, 
which  is  in  agreement  with  Sadik  (1979)  who  used  the  same 
comparison  and  check  stars.  Table  1  summarizes  the 
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catalog  data  on  UV  Piscium,  and  the  comparison  and  check 
stars.  Spectral  types  and  magnitude  information  are 
from  the  H.  D.  catalog  for  comparison  purposes. 


Table  1 
Catalog  Data 


UV  Piscium 

Comparison 

Check 

HD 

7700 

7997 

7918 

BD 

+6°  189 

+6°  197 

+6°  195 

GC 

— - 

1596 

1585 

R.  A. 

(1900) 

lh  11 *1*7 

lh  14?6 

lh  13*1*7 

Dec . 

(1900) 

+6°  16' 

+6°  58' 

+6°  54' 

Ptm 

91?! 

8*1*7 

7*1*9 

Ptg 

9*1*9 

9*?5 

8*1*7 

Spectral  Type 

G5 

G5 

G5 

The  observational  program  was  planned  in  order 
to  obtain  complete  coverage  of  the  light  curve  during 
each  month  of  the  observing  season.  However,  due  to 
limitations  on  the  number  of  scheduled  nights  available 
and  cancelled  nights  because  of  inclement  weather,  total 
monthly  coverage  was  not  achieved,  but  a  single  light 
curve  giving  complete  coverage  of  all  phases  was  obtained. 
The  K-line  filters  were  discontinued  after  the  first 
two  months  in  order  to  increase  the  number  of  individual 
uvby  filter  readings.  This  was  necessitated  mainly 
because  of  the  long  integration  times  required  for  each 


filter  which  severely  limited  the  number  of  individual 
filter  readings  obtainable  each  night. 

Four-color  uvby  standard  stars  were  observed 
on  two  nights,  one  at  the  start  of  the  observing  season, 

6  September,  1980  U.T. ,  and  one  near  the  end,  11  December, 
1980  U.T.  On  the  first  night,  five  standard  stars  were 
observed  repeatedly  throughout  the  night  in  the  order 
standard,  sky,  standard  followed  by  a  comparison,  sky 
sequence  which  was  then  repeated  for  the  next  standard. 

The  filters  were  observed  in  the  sequence  ybvu  with  one 
observation  in  each  filter.  No  averaging  of  the  obser¬ 
vations  was  performed,  and  each  observation  was  treated 
separately.  On  the  second  night  of  standard  observation, 
13  standard  stars  were  used  and  observed  in  the  sequence 
sky,  standard,  sky  followed  by  sky,  comparison,  sky 
sequences.  The  order  of  filter  observation  for  the  stan¬ 
dards  was  yybbwuu  and  the  two  deflections  in  each  filter 
were  averaged  during  reduction  of  the  data. 

Reduction  of  the  observational  data  was  performed 
on  the  San  Diego  State  University  Computer  Center's  IBM 
1130  computer.  The  punched  paper  tape  from  each  night's 
observations  was  converted  to  computer  cards.  These 
were  then  used  as  the  input  to  the  Department  of  Astro¬ 
nomy's  differential  photometry  reduction  program  FOTOM. 

For  each  observation  of  the  variable  star  in  each  filter, 
the  program  computed  the  Heliocentric  Julian  Date,  the 
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magnitude  difference  in  the  sense  of  variable  star  minus 
comparison  star,  and  the  orbital  phase.  The  data  con¬ 
sisted  of  188  observations  in  the  y  filter,  187  in  the 
b,  188  in  the  v  and  188  in  the  u.  In  addition,  83  obser¬ 
vations  for  the  variable  were  obtained  in  both  the  broad 
and  narrow  K-line  filters.  The  K-line  observations  were 
not  used  to  form  delta  magnitudes  but  were  used  to  form 
an  index  in  the  sense  of  narrow  filter  minus  broad  filter. 
A  similar  index  was  formed  for  observations  taken  of 
the  comparison  star. 

The  data  for  the  four  uvby  filters  are  tabulated 
in  Appendix  A  and  includes  the  Heliocentric  Julian  Date, 
the  phase  and  the  magnitude  difference.  The  phase  for 
each  observation  was  calculated  from  the  equation, 

PHASE  =  (H.J.D.-(Eq+EP) )/P,  where  H.J.D.  is  the  Helio¬ 
centric  Julian  Date  of  the  observation,  P  is  the  period, 
and  E  is  the  number  of  cycles  which  have  elapsed  from 
some  base  epoch,  Eq.  The  ephemeris  given  by  Huth 
(1959)  was  used  with  EQ  =  2428038.555,  and  P  =  0.861046 
days . 

FOTOM  was  also  used  to  calculate  color  index  dif¬ 
ferences  of  Mb-y)  and  A(u-b),  again  in  the  sense  of  (b-y) 
for  the  variable  star  minus  the  (b-y)  value  for  the  com¬ 
parison  star.  Both  of  the  indexes  were  calculated  directly 
by  subtracting  the  extinction  corrected  value  of  y  from 
that  of  b  and  similarly  for  (u-b).  It  should  be  noted 
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that  these  are  not  transformed  color  indexes  on  the  stan¬ 
dard  system.  After  completion  of  the  observing  runs, 
the  author  decided,  not  to  transform  the  raw  magnitudes 
to  the  standard  system  because  the  large  integration 
times  required  often  resulted  in  up  to  a  five  minute 
period  between  the  start  of  a  y  filter  observation  and 
the  ending  of  a  u  filter  observation  for  one  sequence 
of  the  variable.  This  fact  in  conjunction  with  the  10 
to  15  minutes  or  more  required  between  each  variable 
star  sequence  would  have  required  extensive  interpolation 
in  order  to  calculate  the  appropriate  standard  indices. 
Such  interpolation  when  considered  with  the  short  period 
of  this  binary  could  begin  to  cover  significant  portions 
of  the  orbital  phase,  particularly  in  the  phase  portions 
during  primary  and  secondary  eclipse.  Any  values  of 
the  indices  obtained  during  these  phase  portions  would 
be  of  questionable  value  at  best.  Even  the  color  indices 
formed  directly  from  the  filter  magnitudes  should  be 
used  with  this  stipulation  that  the  values  may  be  erro¬ 
neous  in  the  areas  of  the  primary  and  secondary  eclipse. 

Since  observations  of  uvby  standard  stars  had 
been  made  previous  to  this  decision,  the  author  still 
elected  to  transform  the  comparison  star  to  the  standard 
system.  The  input  for  the  transformation  was  the  raw 
magnitudes  of  each  filter  and  the  mean  air  mass  for  each 
set  of  individual  observations.  These  values  were 


computed  by  the  FOTOM  reduction  programs  for  each  night's 
observations . 

A  Fortran  program  was  then  written  based  on  the 
method  suggested  by  Crawford  and  Barnes  (1970).  This 
program  was  run  on  the  San  Diego  State  University  Computer 
Center's  IBM  360  computer.  After  the  program  calculated 
the  observed  indices,  (b-y)',  m^ ' ,  and  c^',  for  the  uvby 
system  from  the  raw  magnitudes,  a  set  of  initial  extinc¬ 
tion  coefficients,  K,  ^ »  and  Kj,  was  used  to  calcu¬ 
late  the  extinction  corrected  values  using  the  equations 

yobs  -  y  -  KX 

(b-y}obs  =  (b-y}'  -  v 

(obs)  =  m^ '  -  ^X 

c1(obs)  =  c^'  -  K-jX 

where  X  is  the  mean  air  mass.  These  values  were  then 
used  along  with  the  respective  values  on  the  standard 
system  for  that  star  in  a  least  square  solution  to  deter¬ 
mine  the  transformation  coefficients  from  the  equations 

V  “  A  +  B(b-y)  +  yQbs 

(b-y)  -  C  +  D(b-y)obs 

m^  =  E  +  Fm^(obs)  +  J(b-y) 
ct  =  G  +  Hc^(obs)  +  I (b-y) 

New  extinction  coefficients  were  then  determined  using 
these  transformation  coefficients,  and  the  calculations 
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were  repeated  until  convergence  of  all  of  the  coefficients 
was  obtained. 

An  observational  error  was  inadvertently  made 
on  the  night  of  11  December  1980.  The  standards  were 
all  observed  at  or  near  the  zenith  and  did  not  cover 
a  very  large  range  of  air  mass.  As  a  result,  the  extinc¬ 
tion  coefficients  were  not  well  determined  and,  hence, 
accurate  transformation  coefficients  could  not  be  calcu¬ 
lated. 


The  values  obtained  for  the  transformation  coeffi¬ 
cients  for  the  night  of  6  September  1980,  are  listed 
in  Table  2  along  with  the  mean  values  obtained  for  the 
comparison  star  based  on  these  coefficients. 

Table  2 


Comparison  Star  Transformation  Values 


A 

-1.6412 

I 

-0.3587 

B 

-0.0824 

K 

0.2432 

C 

0.9336 

K1 

0.0790 

D 

1.1065 

k2 

0.0789 

E 

-0.5807 

k3 

0.1854 

F 

1.0789 

V  (S.D. ) 

8.49  (0.06) 

J 

-0.1717 

(b-y)  (S.D.) 

0.61  (0.06) 

G 

0.3490 

(S.D. ) 

0.45  (0.08) 

H 

0.8135 

c^  (S.D.) 

0.68  (0.04) 

Chapter  3 


EPHEMERIS  CORRECTION 

Two  primary  eclipses  were  observed  in  their 
entirety  on  two  of  the  observing  nights  for  UV  Piscium. 

The  phase  of  each  observation  was  obtained  from  the 
ephemeris  given  by  Huth  (1959)  as  was  discussed  in  the 
last  chapter.  When  the  magnitude  difference  versus  phase 
was  plotted  as  in  Figure  3,  it  was  immediately  apparent 
that  the  times  of  primary  and  secondary  minimum  were 
shifted  to  the  right  and  were  not  at  the  phase  points 
0.0  and  0.5  respectively  as  they  should  have  been.  Such 
a  shift  can  normally  be  attributed  to  one  of  two  causes, 
either  a  period  change  or  an  incorrect  ephemeris. 

The  ephemeris  for  an  eclipsing  binary  is  an 
equation  of  the  form 

H. J.D.  (Pri.  Min. )  =  Eq  +  EP 
which  predicts  the  Heliocentric  Julian  Date,  H.J.D.,  of 
the  primary  minimum  from  a  base  date,  EQ,  given  the  number 
of  cycles,  E,  elapsed  from  EQ  and  the  period,  P,  of  the 
orbit.  If  the  ephemeris  is  correct,  the  phase  of  the 
times  of  primary  minimum  should  fall  at  0.0. 

The  ephemeris  of  Huth  (1959) 

H. J.D.  (Pri.  Min.)  =  2428038.555  +  0.861046E 
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was  determined  from  photographic  determinations  of  eclipse 
timings  which  do  not  contain  the  greater  accuracy  of 
today's  photoelectric  observations.  Using  the  ephemeris, 
one  can  calculate  a  predicted  time  of  minimum,  C,  and 
with  the  corresponding  observed  time  of  minimum,  0,  can 
obtain  the  difference,  (0-C),  which  can  be  plotted  against 
either  E  or  H.J.D.  Such  a  plot  can  provide  an  indication 
of  any  changes  which  have  occurred  in  the  period  of  an 
eclipsing  binary  or  the  accuracy  of  the  determined  period. 
If  the  ephemeris  is  accurate  with  no  period  changes  having 
occurred,  the  (0-C)  plots  should  all  fall  along  a  hori¬ 
zontal  line  on  the  graph.  If  a  period  change  occurs, 
it  will  appear  as  a  sudden  change  in  the  slope  of  the 
(0-C)  trend  line,  upward  sloping  for  a  period  increase 
and  downward  sloping  for  a  decrease.  If  the  change  in 
period  is  gradual,  it  will  appear  as  a  continually 
changing  slope  either  up  or  down. 

If,  however,  no  period  change  has  occurred,  but 
the  period  was  incorrectly  determined  to  begin  with, 
then  the  trend  line  will  have  a  constant  slope  across 
the  graph  either  upward  or  downward.  An  upward  sloping 
trend  line  would  indicate  the  original  period  was  calcu¬ 
lated  too  short.  This  appears  to  be  the  case  for  UV 
Piscium. 

A  search  of  the  literature  available  through 
the  San  Diego  State  University  Library  and  its  interlibrary 
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loan  service  yielded  83  timings  of  the  primary  minima  plus 
the  two  from  this  current  observational  program.  Seven 
of  the  photographic  timings  had  (0-C)  values  as  large 
as  one-half  day  in  magnitude.  Such  large  values,  which 
were  on  the  order  of  half  the  orbital  period,  were  con¬ 
sidered  as  being  due  to  possible  observational  errors 
and,  hence,  were  not  included  in  subsequent  calculations. 
The  remaining  timings  are  plotted  as  (O-C)'s  versus  H.J.D. 
in  Figure  1  based  on  the  Huth  ephemeris.  The  photographic 
and  visual  timings  along  with  their  (O-C)'s  and  Huth 
epoch  numbers  are  listed  in  Appendix  B.  The  photometric 
values  are  listed  in  Table  3. 


Table  3 


Photometric  Minima  of  UV  Piscium 


H.J.D. 

(2400000+) 

(0-C) 

Huth 

(0-C) 

This  Study 

Reference 

39388.87432 

+0.00984 

-0.00070 

Carr 

39406.95616 

+0.01082 

+0.00024 

Carr 

39407.81725 

+0.01087 

+0.00029 

Carr 

40156.92849 

+0.01209 

+0.00002 

Oliver 

40466.90642 

+0.01346 

+0.00077 

Oliver 

43053.49255 

+0.01740 

-0.00042 

Asian 

43400.49515 

+0.01846 

-0.00004 

Asian 

43406.52250 

+0.01849 

-0.00003 

Asian 

43425.46595 

+0.01893 

+0.00037 

Asian 

43463.34911 

+0.01607 

-0.00257 

Sadik 

43785.38345 

+0.01920 

-0.00007 

Aslan 

44519.85917 

+0.02268 

+0.00195 

This  study 

44526.74563 

+0.02078 

+0.00003 

This  study 
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The  timings  of  primary  minimum  for  this  study 
as  well  as  those  of  Sadik  (1979)  and  Oliver  (1974)  were 
evaluated  using  a  program  originally  written  by  Brian 
Johnson  for  use  on  the  IBM  1130  computer.  This  program 
solved  for  the  minimum  by  reflection  of  the  light  curve 
about  an  initial  estimate  of  the  time  of  minimum  using 
the  method  of  least  squares.  The  program  was  modified 
slightly  by  this  author,  and  the  values  obtained  for 
the  minima  observed  by  Sadik  were  checked  against  his 
values  and  agreed  to  within  0.0003  days. 

An  inspection  of  the  plot  reveals  a  large  scatter 
in  the  photographic  and  visual  (0-C)  values  especially 
when  compared  to  the  13  photometric  values  which  are 
those  points  after  H.J.D.  2439000.0.  There  are  no  abrupt 
or  gradual  changes  in  period  apparent  from  the  plot, 
but  there  is  a  trend  from  lower  left  to  upper  right  which 
is  readily  discernible  from  the  photometric  data  points. 
The  Huth  ephemeris  was  determined  using  observations 
between  H.J.D.  2427000  and  2435000.  Due  to  the  large 
scatter  in  this  range,  it  is  apparent  that  the  period 
could  have  easily  been  calculated  in  error,  and  it  is 
only  when  the  very  early  observations  and  the  more  recent 
ones  are  included  that  the  difference  in  periods  becomes 
apparent.  Aslan  (1978)  investigated  this  period  differ¬ 
ence  using  photometric  timings  obtained  in  1976,  1977 
and  1978  along  with  the  timings  of  Huth.  He  failed  to 
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obtain  the  very  early  timings  and  attributed  the  differ¬ 
ences  between  the  Huth  data  and  the  recent  observations 
to  a  period  change.  Again  it  is  only  in  the  presence 
of  the  very  early  observations  that  one  realizes  that 
the  difference  lies  in  the  calculation  of  the  original 
period.  Carr  (1967)  attempted  to  rectify  this  difference, 
but  he  only  had  his  three  photometric  timings  available 
plus  six  photographic  and  visual  timings  from  Strohmeier 
and  Krigge  (1960).  He  calculated  an  ephemeris  using 
this  data  of  H.J.D.  =  2433155.754  +  0.861047E.  Aslan 
determined  that  based  on  his  photometric  observations 
an  appropriate  ephemeris  would  be  H.J.D.  =  2428038.034 
+  0.8610482E.  It  should  be  noted  that  both  of  these 
periods  are  longer  than  that  of  Huth  as  one  would  expect 
from  the  (0-C)  diagram. 

Based  on  the  more  complete  data  available,  several 
new  ephemerides  were  calculated  by  a  least  square  fit 
to  the  equation  H.J.D.  =  Eq  +  PE  using  the  epoch  numbers, 
E,  of  Huth.  The  first  ephemeris,  calculated  using  only 
the  photographic  and  visual  data,  was 

H.J.D.  =  2428038.55474  +  0.8610474E 
One  should  note  the  very  slight  difference  in  the  sixth 
decimal  place  of  the  period  between  this  and  the  Huth 
period.  If  Huth  had  had  the  earlier  timings  available, 
his  ephemeris  would  have  likely  been  closer  to  this  value. 

A  second  ephemeris  was  calculated  using  all  of 
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the  timings  available  which  gave 

H.J.D.  =  2428038.55354  +  0.8610472E 
Finally,  a  third  ephemeris  was  calculated  based  only 
on  the  photoelectric  data  which  was 

H.J.D.  =  2428038.54259  +  0.8610477E 
All  three  of  the  periods  differ  only  in  the  seventh  deci¬ 
mal  place,  and  they  all  fall  between  the  two  periods 
determined  by  Carr  and  Aslan. 

A  new  ephemeris  was  then  calculated  using  only 
the  13  photometric  times  of  minimum  and  assigning  new 
epoch  numbers  counted  from  a  base  date  which  was  bracketed 
by  the  photometric  observations.  The  new  ephemeris  is 
H.J.D.  »  2441815.30636  +  0.8610477E 
The  epoch  numbers  and  the  (0-C)  values  from  this  equation 
are  listed  in  Table  3,  and  the  (O-C)'s  are  plotted  against 
H.J.D.  in  Figure  2. 


of  Now  Photometric  liphcmeris 


Chapter  4 


LIGHT  CURVES  AND  THEIR  SOLUTION 


Four  separate  and  complete  light  curves  were 
obtained  for  UV  Piscium,  one  in  each  of  the  ubvy  filters. 
A  solution  for  each  light  curve  was  obtained  using  a 
synthetic  light  curve  solution  method  from  which  an 
adopted  solution  was  derived  for  the  system. 


Light  Curves 

Each  of  the  light  curves  was  plotted  on  a  large 
scale  as  magnitude  difference  versus  phase.  These  plots 
are  reproduced  in  Figures  3,  4,  5,  and  6  on  a  smaller 
scale,  one  for  each  filter.  A  number  of  initial  conclu¬ 
sions  can  be  drawn  from  an  inspection  of  the  light  curves. 

The  orbit  of  the  system  is  very  close  to  being 
circular,  which  can  be  seen  from  the  phase  displacement 
of  the  secondary  minimum  relative  to  the  primary,  just 
one-half  of  the  total  cycle.  If  the  phase  error  was 
removed  as  was  discussed  in  Chapter  3,  the  two  minima 
would  occur  at  phases  of  0.5  and  0.0  respectively. 

The  shape  of  the  eclipse  bottoms  can  reveal  the 


type  of  eclipse  taking  place.  The  primary  is  deep  and 
sharp  indicating  it  is  due  to  either  a  partial  eclipse 
or  a  transit.  It  is  also  much  deeper  than  the  secondary 


Phase 


aouaaajjTQ  ap 


eclipse,  indicating  the  primary  star  is  probably  much 
hotter  than  the  secondary.  The  primary  eclipse  also 
becomes  progressively  deeper  going  from  the  y  filter 
toward  shorter  wavelengths  in  the  u  filter.  The  depths 
range  from  about  0.86  magnitude  to  about  1.1.  Just  the 
opposite  effect  is  observed  in  the  secondary  eclipse, 
again  indicating  the  large  temperature  difference.  Here 
the  depths  range  from  about  0.25  magnitude  in  the  y  filter 
to  about  0.17  in  the  u  filter.  These  depths  of  the 
eclipses  are  in  excellent  agreement  with  those  obtained 
by  Sadik,  Oliver  and  Carr.  The  shape  of  the  secondary 
eclipse  is  not  as  sharply  defined  at  minimum  as  the  pri¬ 
mary  is.  The  scarcity  of  observations  in  each  curve 
at  the  minimum  point  makes  it  difficult  to  definitely 
conclude  if  the  eclipse  is  a  partial  or  complete  eclipse. 
However,  there  does  seem  to  be  an  indication  of  a  flat 
bottom  in  the  u  and  v  filters.  Such  a  flat  bottom  would 
indicate  an  occultation  which  would  suggest,  along  with 
the  inspection  of  the  primary  eclipse,  that  the  primary 
star  is  the  larger  one  in  the  system.  This  would  also 
indicate  that  the  inclination  of  the  orbit  is  very  close 
to  90  degrees.  The  duration  of  an  eclipse  is  about  two 
and  one-half  hours. 

The  eclipses  appear  to  be  well  defined  and  symmet¬ 
rical  about  the  minimum  except  for  the  secondary  in  the 
u  filter,  which  can  be  attributed  mostly  to  scatter  in 
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the  observations.  Both  the  ingress  and  egress  from  the 
eclipse  are  relatively  well  defined.  This  would  indicate 
any  oblateness  in  the  two  stars  is  relatively  small. 
Hence,  there  is  little  tidal  distortion  present  even 
with  the  small  orbit  and  close  proximity  of  the  stars. 
Except  for  a  series  of  observations  in  the  phase  ranges 
from  0.3  to  0.5  and  0.7  to  0.9,  which  will  be  discussed 
later,  both  maxima  are  well  defined,  relatively  flat 
and  of  the  same  magnitude  in  each.  It  appears  that  the 
type  of  assymmetry  noted  by  Carr  and  Sadik  is  not  present 
in  these  observations.  The  observational  scatter  is 
on  the  order  of  0. 3-0.4  magnitude  but  increases  to  about 
0.6-0. 8  magnitude  in  the  u  filter.  Any  reflection  effect 
present  is  relatively  small. 

Figures  7  and  8  show  the  color  curves  A(b-y) 
and  A(u-b)  obtained  from  these  observations.  The  A(b-y) 
curve  shows  almost  no  difference  in  or  out  of  eclipse. 
There  is  a  slight  reddening  of  the  color  during  the  pri¬ 
mary  eclipse  as  would  be  expected,  but  it  is  not  pro¬ 
nounced.  The  Mu-b)  color  curve  shows  a  pronounced 
difference  during  eclipse,  and  provides  a  definite 
indication  of  the  spectral  difference  in  the  two  stars, 
with  the  secondary  being  the  cooler  star. 

Figures  9  and  10  are  the  Calcium  II  K-line  plots 


for  UV  Piscium  and  the  comparison  star  respectively. 
As  was  noted  in  Chapter  2,  the  information  obtainable 


gure  9.  Calcium  11  K-Line  for  UV  Piscium 
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from  the  K-line  is  questionable  due  to  filter  irregu¬ 
larities.  This  is  born  out  by  the  two  diagrams.  There 
is  a  very  large  amount  of  scatter,  and  no  trends  relative 
to  phase  can  be  discerned.  In  addition,  the  data  points 
at  phase  0.4  and  0.7-0. 8,  which  fall  above  the  general 
level  of  the  star  points,  correspond  to  the  similar  sets 
of  points  on  the  comparison  star  plots.  Both  respective 
phase  sets  were  observed  on  the  same  nights.  The  compari¬ 
son  star  phase  values  were  assigned  by  time  interpolation 
between  the  respective  phases  of  variable  observations 
on  each  night.  Hence,  these  variations  can  be  attributed 
to  possible  atmospheric  transmission  differences  from 
night  to  night.  A  plot  of  the  points  by  month  of  obser¬ 
vation  also  revealed  no  apparent  differences. 

The  most  notable  aspect  of  the  light  curves  is 
the  absence  of  the  general  asymmetry  between  the  maxima, 
which  was  noted  previously  by  Carr  and  Sadik.  However, 
there  is  a  complication  present  in  each  light  curve  in 
the  maxima  between  phases  0.30  to  0.45  and  0.70  to  0.80. 
The  scatter  in  the  individual  points  is  greater  in  these 
areas  than  the  rest  of  the  maxima.  Plots  of  the  data 
by  each  month  indicated  that  this  scatter  is  due  to  ver¬ 
tical  shifts  of  the  general  level  of  the  maxima  from 
one  month  to  another.  The  phase  region  0.30  to  0.45 
was  observed  on  three  separate  nights,  each  in  a  different 
month,  14  September,  9  October  and  4  November,  1980. 
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The  first  two  nights'  observations  were  at  the  general 
level  of  the  maxima,  but  the  observations  from  4  November 
were  shifted  vertically  upward  above  the  average  maxima 
level  in  all  four  light  curves.  Similarly,  the  phase 
region  from  0.70  to  0.80  was  observed  on  the  nights  of 
11  September,  6  November  and  14  December,  1980.  Again 
the  observations  from  November  and  December  were  at  the 
general  level  of  the  maxima,  but  the  September  obser¬ 
vations  were  shifted  vertically  above  the  maxima  level 
in  all  four  light  curves.  Carr  also  alluded  to  such 
vertical  shifts  in  his  data  when  plots  for  individual 
nights  were  made.  However,  both  his  and  Sadik's  obser¬ 
vations  were  made  over  a  short  period  of  time,  less  than 
one  month,  and  could  not  show  any  shift  in  the  distortions 
from  month  to  month. 

From  the  A(b-y)  values,  a  mean  level  for  the 
two  maxima  was  determined.  They  were  -0.119  and  -0.120 
magnitude  for  the  first  and  second  maxima  respectively. 

The  standard  deviation  for  both  was  0.014  magnitude. 

These  values  excluded  the  two  nights  of  observations 
where  the  distortions  were  present.  The  mean  Mb-y) 
level  was  then  determined  for  the  two  nights  of  shifted 
observations.  They  were  -0.126  and  -0.123  magnitude 
for  the  first  and  second  maxima  respectively.  Although 
both  of  these  means  fall  within  one  standard  deviation 
of  the  mean  maximum  levels,  the  fact  that  all  of  the 
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observations  of  these  two  nights  are  above  the  general 
level  is  of  significant  interest.  The  values  indicate 
that  the  distortions  are  probably  due  to  an  enhancement 
of  the  light  levels  received  from  the  binary.  Such  an 
enhancement  could  be  caused  by  either  an  increase  in 
brightness  due  to  a  hot  spot  or  to  a  localized  increase 
in  brightness  due  to  a  reduction  in  the  surface  coverage 
by  dark  spots  on  one  or  both  stars. 

The  lack  of  complete  coverage  of  the  light  curve 
from  month  to  month  precludes  any  definite  conclusion 
on  movement  of  the  distortion  in  phase  along  the  curve. 
However,  an  inspection  of  the  monthly  observations  of 
different  parts  of  the  light  curve  relative  to  the  posi¬ 
tion  of  the  distortions  suggested  a  retrograde  movement 
in  phase  of  the  distortion  along  the  light  curve.  Such 
a  movement  toward  decreasing  phase  is  a  characteristic 
of  the  wave  distortions  in  RS  Canum  Venaticorum  type 
light  curves  but  on  a  much  shorter  time  scale. 

Solution  of  the  Light  Curve 

The  light  curve  solutions  were  obtained  using 
a  synthetic  solution  method.  The  method  is  based  on 
the  model  developed  by  Nelson  and  Davis  (1972).  Correc¬ 
tions  and  additions  to  that  model  were  made  by  Etzel 
(1975)  and  incorporated  into  a  computer  program  EBOP 
(Eclipsing  Binary  Orbit  Program).  Further  revisions 


of  the  model,  now  referred  to  as  the  Nelson-Davis-Etzel 
model  (NDE) ,  are  contained  in  version  12  of  the  program 
which  was  used  in  this  solution.  The  program  was  run 
on  the  IBM  360  computer. 

The  NDE  model  has  been  demonstrated  by  Popper 
and  Etzel  (1981)  as  an  excellent  model  for  the  analysis 
of  the  light  curves  of  well-detached  eclipsing  binaries. 

UV  Piscium  is  such  a  system  although  some  small  oblateness 
of  the  components  due  to  their  close  proximity  can  be 
expected.  The  NDE  model  incorporates  the  effects  of 
oblateness  after  that  of  Binnendijk  (1960)  for  a  simple 
biaxial  (rotational)  ellipsoid  and  adds  these  effects 
into  a  simple  spherical-limb-darkened  model  to  match 
the  observations.  Reflection  is  treated  essentially 
the  same  as  in  the  simple,  uniformly  illuminated  hemi¬ 
sphere  model  given  by  Binnendijk  (1960).  The  amount 
of  reflection  effect  is  also  assumed  to  be  small  in  the 
NDE  model. 

The  model  parameters  used  in  the  solution  are 
listed  in  Table  4.  Initial  values  for  the  parameters 
were  selected  from  those  of  Sadik  (1979).  The  value 
of  the  radius  of  the  primary  star,  r  ,  is  in  terms  of 
the  semi-major  axis  of  the  orbit.  Theoretical  values 
of  limb  darkening  were  interpolated  from  the  tables  by 
Grygar,  et  al.,  (1972)  for  each  filter,  assuming  main 
sequence  G2  and  K0  stars.  Gravity  darkening  was  set 
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1 


Symbol 


ecosui 


esinw 


Table  4 

Model  Parameters 


Definition 


Central  Surface  Brightness  of  Secondary  (Jp  = 

Equivalent  Spherical  Radius  of  Primary 

Ratio  of  Radii  (spherical,  k  »  r  /r  ) 

s  p 

Limb  Darkening  of  Primary 
Limb  Darkening  of  Secondary 
Orbital  Inclination  (Degrees) 


Orbital  Eccentricity  and  Longitude  of  Periastron 

Gravity  Darkening  of  Primary 

Gravity  Darkening  of  Secondary 

Reflected  Light  from  Primary 

Reflected  Light  from  Secondary 

Mass  Ratio  (q  =  M  /M  ) 

s  p 

Tidal  Lead/Lag  Angle  of  Stellar  Axes 

Third  (Extra)  Light 

Phase  Correction  for  Epoch  Error 


5R£f>- 


SFACT  Luminosity  Scaling  Factor 
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at  zero  due  to  the  small  oblateness  of  this  system  and 
late-type  stars.  Unit  weights  were  used  for  all  obser¬ 
vations.  The  value  of  q  used,  0.75,  is  that  given  by 
Popper  (1976).  Both  t  and  were  set  to  zero.  The 
phase  correction  as  was  calculated  for  each  light  curve 
to  shift  the  primary  minimum  to  phase  0.0.  The  variable 
SFACT  is  the  luminosity  normalization  factor  which  equates 
the  brightness  of  the  theoretical  light  curve  to  the 
observed  curve. 

The  adopted  solution  in  each  filter  was  arrived 
at  after  a  large  number  of  computer  runs  in  each  filter. 
The  primary  criteria  for  solution  acceptance  were  conver¬ 
gence  of  any  one  solution  and  minimization  of  the  sum 

2 

of  the  square  of  the  residuals,  L(o-c)  .  Initial  solu¬ 
tions  were  obtained  for  ecosw  and  esinw.  The  values 
of  e  obtained  indicated  e  was  less  than  0.01.  Hence, 
all  subsequent  solutions  were  made  with  ecosw  and  esinw 
set  to  zero. 

Solutions  were  then  obtained  using  various  com¬ 
binations  of  the  parameters  as  variables  while  others 
were  held  constant  until  the  adopted  solutions  were 
obtained  in  each  filter.  These  values  are  listed  in 
Table  5  along  with  their  standard  errors  when  available. 
Also  listed  are  the  oblateness  e,  semi-major  axis  a  and 
semi-minor  axis  b  of  each  star,  the  fractional  system 
luminosity  L  of  each  component  at  quadrature,  and  the 
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Solution 

Table  5 

Values  and  Standard  Errors 

Variable 

y 

b 

V 

u 

Js 

0.383 

.012 

0.357 

.014 

0.262 

.015 

0.221 

.021 

r 

0.2422 

0.2475 

0.2519 

0.2427 

P 

.0030 

.0032 

.0036 

.0055 

rs 

0.1866 

0.1877 

0.1901 

0.1867 

k 

0.7707 

.0080 

0.7585 

.0146 

0.7547 

.0085 

0.7691 

.0138 

UP 

0.61 

0.70 

0.82 

0.86 

us 

0.72 

0.85 

0.96 

1.00 

i 

88.13 

.50 

89.32 

1.56 

89.15 

1.34 

88.45 

1.00 

SP 

0.0041 

0.0039 

0.0030 

0.0025 

Ss 

0.0112 

0.0116 

0.0124 

0.0122 

AG 

-0.02307 

.00034 

-0.02339 

.00037 

-0.02353 

.00041 

-0.02346 

.00058 

e 

P 

0.0158 

0.0168 

0.0177 

0.0159 

e 

S 

0.0129 

0.0131 

0.0136 

0.0129 

3P 

0.2448 

0.2503 

0.2549 

0.2453 

as 

0.1883 

0.1894 

0.1919 

0.1883 

bp 

0.2409 

0.2461 

0.2504 

0.2414 

bs 

0.1858 

0.1869 

0.1893 

0.1859 

LP 

0.8218 

0.8392 

0.8776 

0.8911 

L. 

0.1782 

0.1608 

0.1224 

0.1089 

S.E.  1 


0.0147 


0.0156 


0.0185 


0.0275 
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standard  error  of  one  observation  s.e.l.  Table  6  lists 
the  mean  values  obtained  from  the  adopted  solutions. 
These  values  agree  closely  to  those  obtained  by  both 
Carr  and  Sadik. 


Table  6 


Mean  Values 

and  Standard 

Deviations 

Parameter 

Value 

S.D. 

ep 

0.0166 

0.0009 

ap 

0.2488 

0.0048 

bp 

0.2447 

0.0045 

E  s 

0.0131 

0.0003 

as 

0.1895 

0.0017 

bs 

0.1870 

0.0016 

i 

88.76 

0.57 

rp 

0.2461 

0.0046 

rs 

0.1878 

0.0016 

k 

0.7633 

0.0079 

Two  further  checks  were  made  on  the  adopted 
solution.  One  was  a  check  on  the  limb  darkening  coeffi¬ 
cients.  The  solutions  were  tested  for  incremental  changes 
in  limb  darkening  away  from  the  theoretical  values. 

No  significant  change  occurred  in  the  residuals  as  mea- 

o 

sured  by  the  £(o-c)  's;  hence,  it  was  concluded  that 
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the  solutions  were  not  sensitive  to  variations  in  limb 
darkening. 

A  test  was  also  made  to  check  the  solution  against 
changes  in  k.  The  value  of  k  was  ranged  from  0.70  to 
0.82  in  order  to  check  the  effect  on  the  residuals. 

In  the  range  from  0.74  to  0.78  for  k,  the  residuals 
increased  less  than  5%,  and  remained  less  than  l07o 
different  in  the  range  from  0.73  to  0.79.  This  suggests 
that  the  true  value  of  k  is  probably  in  the  range  from 
0.74  to  0.78,  which  is  in  excellent  agreement  with  the 
more  formal  standard  error  of  k  in  Table  5  but  gives 
a  more  definitive  range  within  which  k  is  located. 


Chapter  5 
CONCLUSION 

The  absolute  parameters  of  the  orbit  and  stars 
can  now  be  obtained  from  the  adopted  solution  values 
listed  in  Table  6.  Popper  (1976)  gives  the  minimum  masses 
for  the  system,  M^sin  i,  as  1.2MQ  for  the  primary  and 
0.9Mq  for  the  secondary,  which  he  determined  from  spectro¬ 
scopic  observations.  Using  the  value  of  i  from  Table 
6,  the  absolute  masses  were  calculated  and  are  listed 
in  Table  7. 

Table  7 

Absolute  Parameters 

Parameter  Value 


Mp 

1.201 

Mo 

Ms 

0.901 

M 

o 

a 

4.879 

R 

o 

RP 

1.201 

R 

o 

Rs 

0.916 

R 

o 

Rcrit(p) 

1.971 

R 

o 

Rcrit ( s ) 

1.737 

R 

o 
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Table  7  also  lists  the  values  for  the  radius 

of  the  orbit,  a,  and  the  stellar  radii,  R  and  R  .  The 

p  s 

orbital  radius  can  be  found  from  Kepler's  third  law 

a3  =  (74.5445)P2(M1  +  M2 ) 

where  the  units  are  solar  radii  for  a,  days  for  P  and 
solar  masses  for  and  M2 .  The  individual  stellar  radii, 
R y  in  terms  of  solar  radii  can  then  be  obtained  using 
the  relation 


With  these  previously  obtained  values,  it  is 
now  possible  to  determine  the  size  of  the  critical  Roche 
lobes  for  this  system.  Plavec  (1970)  gives  the  following 
algorithm  for  determining  the  Roche  lobes  given  the  mass 
ratio  of  the  system: 


rcrU(j)  *  °'38  +  0-201°SIH^) 


where  rcrit(j)  is  in  units  of  the  separation  of  the  stars, 
a.  The  values  of  the  Roche  lobes  in  absolute  units  are 
given  in  Table  7  and  can  be  compared  to  the  respective 
stellar  radii.  Both  stars  fill  less  than  75%  of  their 
respective  Roche  lobes.  The  primary  and  secondary  fill 
61%  and  53%  of  the  respective  Lobes.  Hence,  this  system 
is  a  well  detached  system. 

There  is  an  ambiguity  as  to  the  spectral-lumi¬ 
nosity  classification  of  this  system,  especially  of  the 
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secondary  star.  The  spectral  class  of  the  primary  has 
been  determined  by  Popper  (1976)  as  G2 ,  and  the  secondary 
was  classified  on  the  basis  of  light  curve  determined 
colors  by  Oliver  as  a  KO.  As  a  check  on  these  values,  an 
estimate  of  the  stellar  colors  was  made  from  this  study, 
but  it  must  be  emphasized  that  these  are  not  definitive 
due  to  problems  referred  to  earlier  in  transforming  the 
colors.  However,  with  this  stipulation  in  mind,  the  (b-y) 
color  on  the  standard  system  of  the  primary,  determined 
from  that  of  the  midpoint  of  the  secondary  eclipse  when 
only  light  from  the  primary  is  observed  is  0.46.  That  of 
the  secondary  star,  although  contaminated  slightly  by  the 
primary  during  the  primary  eclipse,  is  0.52.  These  values, 
when  compared  to  those  given  by  Popper  (1980)  for  main 
sequence  stars,  place  the  primary  as  a  late  G  and  the 
secondary  as  an  early  K  star.  If  the  primary  is  in  fact  a 
G2  star,  it  may  be  slightly  evolved  above  the  main 
sequence,  but  it  definitely  is  not  a  subgiant  as  can  be 
seen  by  its  small  radius.  Similarly,  the  secondary  may  be 
an  early  K  star,  but  it  also  is  definitely  not  a  subgiant 
due  to  its  small  radius  as  was  alluded  to  in  Chapter  1. 
Hence,  this  system  should  be  classified  as  two  main 
sequence  stars  with  the  primary  a  G2  and  the  secondary  an 
early  K  star. 

Two  significant  questions  concerning  the  wave  dis¬ 
tortion  in  the  light  curve  remain  to  be  resolved.  The 
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first  concerns  whether  the  distortion  when  present  is 
really  an  enhancement  of  the  light  output.  This  question 
should  be  addresed  by  an  attempt  to  specifically  define  a 
non-distorted  level  for  the  out  of  eclipse  maxima.  Any 
distortion  then  observed  could  be  referenced  to  a  well  de¬ 
fined  non-distorted  maximum  level  to  determine  if  it 
really  is  an  enhancement. 

The  second  question  concerns  the  movement  of 
the  distortion  from  month  to  month.  Is  this  a  real  move¬ 
ment  of  one  distortion  wave  or  the  appearance  at  different 
phase  positions  of  several  individual  distortions?  Also, 
is  the  motion  prograde  or  retrograde?  Both  of  these 
areas  of  investigation  could  be  addressed  by  an  observing 
program  designed  to  observe  the  entire  light  curve  during 
a  two  week  period  or  less  and  then  comparing  several 
such  curves  obtained  in  one  observing  season,  with  partic¬ 
ular  emphasis  on  detailed  observations  of  the  maxima. 

Such  a  program  should  also  be  carried  out,  if  possible, 
on  a  telescope  larger  than  the  41-cm  telescope  used  in 
this  study  due  to  the  magnitude  range  of  UV  Piscium. 

In  conclusion,  UV  Piscium  has  been  found  to  be  a 
well  detached  binary  composed  of  two  main  sequence  stars 
of  spectral  class  G2  and  early  K.  The  light  curves  are  well 
defined  with  very  little  distortion,  except  for  a  small 
wave  distortion  in  the  maxima,  and  their  solutions  have 
permitted  reliable  absolute  elements  to  be  determined. 
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